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ABSTRACT 

Context. Instrumental projects that will improve the direct optical finding and characterisation of exoplanets have advanced sufficiently 
to trigger organized investigation and development of corresponding signal processing algorithms. The first step is the availability of 
field-of-view (FOV) models. These can then be submitted to various instrumental models, which in turn produce simulated data, 
enabling the testing of processing algorithms. 

Aims. We aim to set the specifications of a physical model for typical FOVs of these instruments. 

Methods. The dynamic in resolution and flux between the various sources present in such a FOV imposes a multiscale, independent 
layer approach. From review of current literature and through extrapolations from currently available data and models, we derive the 
features of each source-type in the field of view likely to pass the instrumental filter at exo-Earth level. 

Results. Stellar limb darkening is shown to cause bias in leakage calibration if unaccounted for. Occurrence of perturbing background 
stars or galaxies in the typical FOV is unlikely. We extract galactic interstellar medium background emissions for current target lists. 
Galactic background can be considered uniform over the FOV, and it should show no significant drift with parallax. Our model 
specifications have been embedded into a Java simulator, soon to be made open-source. We have also designed an associated FITS 
input/output format standard that we present here. 

Key words. Instrumentation: high angular resolution - Methods: analytical - astronomical data bases: miscellaneous - astrometry - 
ISM: structure - Galaxy: stellar content 



1. Introduction 

Instruments currently under design for direct optical exoplane- 
tary search and characterisation need to go beyond the indirect 
techniques used so far for the discovery of the ~ 200 currently 
known exoplanets, and must collect planetary photons. Beyond 
the joint determination of the science objectives of albedo, plan- 
etary radius and orbital parameters, the major aim of these in- 
struments is to establish the presence, in a potential atmosphere, 
of chemical markers of life processes (biomarkers). This would 
be done through the detection of their absorption features in the 
spectral flux emitted by the planet. 

Because of this, broadband observation is required. Two 
types of spaceborne instruments are currently under develop- 
ment. Both types of instruments reject stellar light, so that the 
10~ 9 and 10" 6 respectively weaker planetary flux is detectable in 
the residual noise. The Terrestrial Planet Finder Coronograph 
(TPF-C) is a 6 x 2 m 2 monolithic collector space telescope in the 
visible (Traub et al. 2006). Free-flying-collector interferometers, 
in the infrared (band extending from 6 to 18 fun), used in a par- 
ticular optical design called nulling interferometry (Bracewell 
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1978), are also being considered. In this latter technique, the op- 
tical array is phased so that light from the on-axis star is de- 
structively interfered. As the array is rotated, the off-axis planets 
pass through the peaks and valleys of the instrumental response 
on the sky (the so-called transmission map), which generates a 
modulated signal. The main interferometric projects are Darwin 
(Leger et al. 1996; Fridlund 2000) and TPF-I (Beichman et al. 
1999). Complementarity of biomarkers at these two wavelength 
ranges, associated with the advantages and shortcomings of each 
of these classes of instruments, explain this parallel effort. 

Both approaches are currently mature enough to trigger or- 
ganized investigation and development of signal processing al- 
gorithms for planet detection and characterisation: Ferrari et al. 
(2006) for direct imaging, and Mugnier et al. (2006), Thiebaut 
& Mugnier (2006), Thiebaut et al. (2007), Marsh et al. (2006), 
Draper et al. (2006) for nulling interferometry. 

This paper specifies a physical and mathematical model of 
source FOVs, called Origin. As will be seen in the various sec- 
tions of this paper, there is an abundance of available elements 
characterizing exoplanetary FOVs. We felt there was a need for 
an integrated access to this information for simulation input, data 
exchange and outreach. 

At the present time, the instruments capable of exo-Earth de- 
tection and characterisation are in a very early definition phase: 
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no concepts are considered final. For this reason we make no 
simplifying assumption regarding the instrument, in particular 
its sensitivity and/or its ability to discriminate between specific 
scene features and/or noise sources. 

This work has been greatly inspired by the European Darwin 
mission, hence the nulling interferometry point of view is often 
significantly developed beyond the conclusions that apply more 
generally to exo-Earth finding instruments. 



2. Framework 

In this section we present the framework elements of our model: 
the building-block rationale and considerations on spatial, spec- 
tral and temporal resolutions. 



2.1. Building-Block Model 

Our simulator of astronomical scenes aims at modeling the an- 
gular and spectral distribution of light received from the ob- 
served exoplanetary system. The proposed model is built by su- 
perimposing the emission of the various sources that are seen 
by the instrument. Following this, the specific intensity (units: 
W rrr 2 rad~ 2 \\mT l ) observed in a direction s is: 



/(s, A, = 7*(s,A,0 + yV/s,A,0 



+ 7 ez (s, A, f) + I Z (A, t) + 7 bg (s, A, t) ., 



(1) 



where 7+ is the star's emission (see Sect. 3), Ij is the contribu- 
tion of the /'-th planet (see Sect. 4), 7 ez is the emission by the 
exozodiacal cloud (see Sect. 5.2), 7 Z is the local zodiacal cloud 
emission assumed to be uniform across the field of view (see 
Sect. 5.1) and, finally, 7b g accounts for the contribution of back- 
ground sources (Sect. 5.3 and Sect. 5.4). If (x,y) are the cartesian 
coordinates of a given source in a plane perpendicular to the line 
of sight, then the angular direction of observation with respect to 
the center of the field of view is s ^ {8 X , 8 y ) = (x,y)/d, where d 
is the distance from the observer. 

For each source type, the various sections of this paper will 
present a discussion on the astrophysical features likely to pass 
the instrumental filter, at the level of the signal that an exo- 
Earth would produce, and hence which require modeling; sub- 
sequently, example layer-outputs for that particular source are 
demonstrated. 



2.2. Coordinates Systems 

Our model accounts for the observatory's Solar-system coordi- 
nates (ecliptic, Earth or L2 1 ) and for target coordinates, enabling 
local zodiacal drift computing (see Sect. 5.1). The planed verti- 
cal extent (under 7x 10 5 km ) of the observatory's halo Lissajous 
libration orbit around L2 (Landgraf 2004) is considered small 
compared to the zodiacal cloud's thickness, and was not imple- 
mented. 

A target (i.e. exosystem)-proper coordinate system is imple- 
mented for calculation of the stellar flux received and phase- 
reflected by a planet (Sect. 4.2). It also enables consistent scene 
generation for revisits, astrophysical community databases im- 
port/export, and robustness testing of image reconstruction al- 
gorithms. 



Temporal accounting has two aspects. Time is used to model 
variations in the scene that occur typically on timescales com- 
parable to those of an observation, such as the motion of close, 
short revolution-period planets (Pegasides), or stellar variability 
(flares). The resolution is typically 1/1 th of an hour. Dates are 
used to initialize scene consistently, for simulating revisits; the 
same 1/1 th hour resolution applies, but expressed in fractional 
Julian Date (JD). In the output of our model, we assume that 
the exposure duration Af is sufficiently short to consider that the 
scene is static during a given exposure. 

2.3. Spatial Resolution 

In the output of the proposed building-block model (see 
Appendix B), sources with apparent size larger than the instru- 
mental resolution must be considered as resolved, and modeled 
by maps of their brightness distribution. These maps (or images) 
must be sampled with a pixel size AO smaller than the instrumen- 
tal resolution limit: 

AG <«; - — (2) 

Umax 

where A m \ n is the smallest spectral channel effective wavelength, 
and 7? max the largest interferometric baseline. For a monolithic 
telescope, 7? max is the diameter of the pupil. 



2.4. Spectral Resolution 

To avoid loss of coherence in the interferometer model due to 
the finite spectral bandwidth AA in the scene model, the spectral 
resolution must be chosen so that the phase difference across a 
spectral channel is negligible: 



Be 



B9 



AA A+kAA 



BOAA 



A 2 



« 1 



(3) 



where A is the effective wavelength of the spectral channel, B is 
the interferometric baseline, and is the view angle. Hence the 
spectral resolution must be chosen so that: 



T^max 0] 



max ^max 



A/1 



A 



(4) 



where max is the radius of the field of view. 

The interferometric baseline is typically chosen so as to have 
the dark fringe of the nulling interferometer not larger than the 
inner size 0Hz in of the habitable zone, that is: A/B msix ~ 0uz ia - 
Assuming a typical field of view radius of 6< max ~ 075, and since 
the smallest considered 6nz m is 0701 mas (Kaltenegger et al. 
2006), the spectral resolution must be better than ~ 50. 

2.5. Photon Counts 

To produce an image, it is convenient to first consider the ex- 
pression of the photon count received into a spectral channel by 
a pixel during an exposure. For a resolved source this is: 



Absolved = m t~ 7(s, A, t) AO 2 AA At , 
he 



(5) 



Second Lagrange point of the Earth-Sun 2-bodies system 



where Af is the exposure duration, S te \ is the collecting area [e.g. 
unobstructed surface of the telescope primary mirror) and t]\ is 
the instrumental throughput at the spectral channel wavelength. 
The photon count in Eq. (5) assumes that the spectral bandwidth 
AA, the pixel size AO, and the exposure Af are small enough 
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compared to typical scales of variation of the specific intensity 
I(s, A, t), as discussed in the previous subsections. 

For an unresolved source (e.g. a planet or any point-like 
source), the term I(s, A, t) AO 2 in Eq. (5) must be replaced by the 
specific flux F{A, t). For the planets, the specific flux can be com- 
puted straightforwardly from atmospheric models (Sect. 4.1) 
which are used as input databases in our building-block model. 
The number of photons received by a pixel is thus: 



N, 



unresolved 



J] a <Std -j — F(A, t) AA Af . 
he 



(6) 



The building-block output of our model comprises two cate- 
gories, accordingly. For resolved sources such as the exozodiacal 
dust emission, the output consists of 3-D (9 X , 9 y , A) maps of the 
brightness distribution of the source at a given observing time t, 
and integrated by a pixel of angular size A6 and effective spectral 
bandwidth AA: 

/resolved(S, A, t) = A0 2 M (?) 

tl C 

in units of number of incident photons per m 2 per s. For un- 
resolved sources, such as the planets, the building-block model 
output is simply the specific flux of the point-like source inte- 
grated in every spectral channel at a given observing time t: 



AF(AJ) 

/unresolved (A, t) — ; AA 



ll C 



in the same units as / r 



resolved- 



3. Star Model 

This section describes the modeling of /^(s, A, ?), the specific in- 
tensity emitted by the star. 



3. 1 . Limb Darkening and Stellar Leakage 

Due to the finite extension of the star, coupled with instrumental 
instabilities, the stellar light may not be completely suppressed 
by the coronagraphic technique, be it classical or interferomet- 
ric coronagraphy. The residual light is called leakage. We would 
like to know if choosing not to model the stellar limb darkening 
would lead to significant errors of the leakage that an instrumen- 
tal simulator would produce. 

Assuming the star has spherical symmetry, the angular and 
spectral distribution of its emission is given by: 



(9) 



where /" = I+Qi = 0, A) is the specific intensity emitted by the 
star in a direction normal to its surface, and D A (p) is the limb- 
darkening law which is a function of the cosine fi of the angle 
between the viewing direction and the normal to the surface: 



ti= Vi - (0/e*) 2 



(10) 



where 6 is the angular direction with respect to the center of 
the star and 0* is the apparent radius of the star. Following 
Van Hamme (1993), we consider the following possible limb- 
darkening laws: 



D A (yu) = 



1 
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-x x {\ 




1 




-») 


1 


-x A (l 


-») 



(black body) 
(linear) 
(logarithmic) 
(square root) 



(11) 



where x A and y A are tabulated parameters which have been com- 
puted by Van Hamme (1993) for 410 stellar spectra synthesized 
with the ATLAS code. 

Now we consider the case of a nulling interferometer. The 
stellar contribution in its output signal reads: 



A+(A,t) 



I 



R(s,A,t)h(s,A,t)ds 



(12) 



where R(s, A, t) is the instrumental instantaneous transmission 
map. At the center of the field of view, it can be approximated 
(Absil 2001), in the ideal, unperturbed case, by a power law. 
We also consider an azimuthal symmetry (similar results are ob- 
tained in a full two-dimensional integration) which leads to: 

R(s,A,t)cc6" . (13) 

We are now going to take this pattern as a model for the ef- 
fective transmission over a realistic exposure, with perturbations. 
Lay (2004, Table 5) shows that, in total, power-law "geometric 
type" leakage contributions are dominant over the "floor type" 
contributions. Indeed, we are not interested in the exact evalua- 
tion of the leakage, but in a useful approximation of the relative 
influence of modeling the integration of the limb-darkened emis- 
sion through a "null profile". 

Then, since 6 = 6+ -y/l - /i 2 and assuming, for sake of sim- 
plicity, a linear limb-darkening law, the star leakage scales as: 



(8) A+(A,t) 



J (\-H 2 f 2 [l-(l-/i)jc a ]j"dp 



1 - y(n) x A 
n + 2 



with 



y(n) = 1 - 



n + 



2 r(I)r(^) 



2 r(¥ ) 

and where F is the gamma function: 



r(z) 







f~ x e~' df . 



(14) 



(15) 



(16) 



Table 1. Numerical values of y(n) defined in Eq. (15) for differ- 
ent nulling power n. 



n: 2 4 6 

y(n) : 0.467 0.543 0.594 



To assume that the star emits as a black body is equivalent to 
taking x A = 0. Thus the term y(n)x A in Eq. (14) is equal to the 
relative attenuation of stellar leakage due to the limb-darkening 
law, compared to a pure black body star. Table 1 shows that 
y(ri) ^ 0.5 for typical values of n, and Table 2 lists the aver- 
age value of x A in the photometric bands. From these tables, we 
can estimate that the relative attenuation of stellar leakage due 
to taking into account the limb-darkening law is y(n) x A ~ 5 % 
in the infrared for F to M stars, but can be as high as ~ 35 % in 
the visible. These figures are consistent with similar analysis by 
Absil et al. (2006). Although computed for a linear law, the same 
results would be obtained for the other limb-darkening laws con- 
sidered and, to summarize, accounting for the limb darkening is 
important so as not to overestimate the stellar leakage in the visi- 
ble range, but could be neglected, in a first approximation, in the 
infrared. Finally it is worth noting that dynamic errors over the 
duration of an exposure will result in more stellar leakage, some 
of which can not be easily calibrated (Lay 2004). 
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Table 2. Mean linear limb-darkening coefficient Xa in photometric bands for stars of different spectral types. 



Type 


U 


B 


V 


R 


I 


J 


H 


K 


L 


M 


N 


Q 


FOV 


0.55 


0.57 


0.48 


0.37 


0.29 


0.23 


0.17 


0.15 


0.13 





11 


0.06 


0.04 


F2V 


0.58 


0.59 


0.49 


0.38 


0.30 


0.24 


0.18 


0.16 


0.13 





11 


0.06 


0.04 


F5V 


0.67 


0.64 


0.52 


0.42 


0.33 


0.27 


0.20 


0.17 


0.14 





12 


0.06 


0.04 


F8V 


0.72 


0.67 


0.55 


0.44 


0.35 


0.28 


0.21 


0.18 


0.15 





12 


0.07 


0.05 


GOV 


0.78 


0.71 


0.57 


0.46 


0.37 


0.30 


0.22 


0.19 


0.16 





13 


0.07 


0.05 


G2V 


0.84 


0.75 


0.61 


0.49 


0.40 


0.32 


0.24 


0.21 


0.17 





13 


0.07 


0.05 


G5V 


0.84 


0.75 


0.61 


0.49 


0.40 


0.32 


0.24 


0.21 


0.17 





13 


0.07 


0.05 


G8V 


0.89 


0.79 


0.64 


0.52 


0.42 


0.34 


0.25 


0.22 


0.17 





14 


0.08 


0.05 


KOV 


0.94 


0.83 


0.68 


0.55 


0.45 


0.36 


0.27 


0.23 


0.18 





14 


0.08 


0.05 


K1V 


0.97 


0.87 


0.73 


0.58 


0.47 


0.38 


0.28 


0.25 


0.19 





14 


0.08 


0.05 


K2V 


0.97 


0.87 


0.73 


0.58 


0.47 


0.38 


0.28 


0.25 


0.19 





14 


0.08 


0.05 


K3V 


1.00 


0.91 


0.76 


0.62 


0.49 


0.41 


0.30 


0.26 


0.20 





15 


0.09 


0.05 


K4V 


1.01 


0.95 


0.80 


0.64 


0.51 


0.43 


0.31 


0.27 


0.21 





16 


0.09 


0.05 


K5V 


0.98 


0.94 


0.79 


0.66 


0.51 


0.43 


0.32 


0.27 


0.21 





16 


0.09 


0.06 


K7V 


0.85 


0.82 


0.70 


0.61 


0.46 


0.37 


0.31 


0.25 


0.18 





14 


0.09 


0.06 


MOV 


0.73 


0.68 


0.61 


0.56 


0.40 


0.29 


0.24 


0.19 


0.14 





13 


0.08 


0.06 


M1V 


0.73 


0.68 


0.61 


0.56 


0.40 


0.29 


0.24 


0.19 


0.14 





13 


0.08 


0.06 


M2V 


0.73 


0.66 


0.62 


0.57 


0.40 


0.25 


0.18 


0.15 


0.12 





12 


0.07 


0.06 


M3V 


0.73 


0.66 


0.62 


0.57 


0.40 


0.25 


0.18 


0.15 


0.12 





12 


0.07 


0.06 


M4V 


0.72 


0.63 


0.58 


0.54 


0.37 


0.23 


0.17 


0.14 


0.11 





11 


0.07 


0.06 



3.2. Spectral Classification 

In the implementation of our scene model, we make use of an 
input database of stellar parameters, for different stellar spectral 
types. This database is derived from the work of Van Hamme 
(1993), who computed the detailed coefficients for various limb- 
darkening laws of synthetic spectra from the ATLAS code. 
These spectra were simulated for solar chemical composition 
stars with a wide range of effective temperatures T+, and sur- 
face gravities g+, covering most of the observed HR diagram. 

By using spectral classification tables (Schmidt-Kaler 1982), 
we derived the MK spectral type and luminosity class of the stars 
from their physical parameters g+ and T+. The classification also 
yields the stars' absolute luminosity L*, and an estimate of their 
mass M+ which is required to compute the orbits of the planets. 
Note that the MK classification only gives average or typical val- 
ues of physical stellar parameters for a given spectral type; these 
values are therefore not consistent with those of an individual 
star. To reduce these inconsistencies, we tuned the parameters so 
that the star radius and surface gravity verify: 



Table 3. Star model parameters in the Origin database for spec- 
tral types F through M. 



and 



R+ — Re 



g* = go 



1/2 /T v2 
* G 



A' 

R* 



(17) 



(18) 



Table 3 lists the resulting star parameters in our input database 
for spectral types F through M. 

3.3. Stellar Features 

Woolf et al. (1998) were the first to devise a technique called 
"chopping", or "modulation", enabling nulling interferometers 
to see only the shot noise of centrosymetric sources in the FOV 
(to a first order stellar leakage, and exozodiacal dust emission). 
Given the 10 6 brightness dynamic between the star and an exo- 
Earth, stellar features may introduce biases in both detection and 
spectroscopy that cannot be removed by means of internal mod- 
ulation. Such problems are expected from the contribution of 
spatially non-symmetric emission features such as non-uniform 



Type 


r* 


MJM e 


LJL 




R*/Ro 


logte*/£o) 


FOV 


7250 


1.60 


5.61 




1.50 


-0.15 


F2V 


7000 


1.52 


4.35 




1.42 


-0.12 


F5V 


6500 


1.40 


2.72 




1.30 


-0.08 


F8V 


6250 


1.18 


1.90 




1.18 


-0.07 


GOV 


6000 


1.05 


1.41 




1.10 


-0.06 


G2V 


5750 


9.96- 10" 1 


1.03 




1.02 


-0.02 


G5V 


5750 


9.20- 10-' 


8.35 ■ 10- 




9.20 ■ 10" 1 


0.04 


G8V 


5500 


8.40- 10-' 


6.35 ■ 10- 




8.77 ■ 10 _1 


0.04 


KOV 


5250 


7.90- 10-' 


4.95 • 10- 




8.50- 10" 1 


0.04 


K1V 


5000 


7.64- 10-' 


3.81 ■ 10" 




8.22 ■ 10" 1 


0.05 


K2V 


5000 


7.40- 10-' 


3.57 ■ 10" 




7.95 ■ 10" 1 


0.07 


K3V 


4750 


7.16- 10-' 


2.72- 10" 




7.69 ■ 10-' 


0.08 


K4V 


4500 


6.92- 10-' 


2.05 • 10- 




7.44 ■ 10" 1 


0.10 


K5V 


4250 


6.70- 10-' 


1.53- 10" 




7.20- 10" 1 


0.11 


K7V 


4000 


6.01 ■ 10-' 


1.03- 10" 




6.69 ■ 10~' 


0.13 


MOV 


3750 


5.10- 10-' 


6.42 ■ 10" 


2 


6.00- 10" 1 


0.15 


M1V 


3750 


4.52- 10-' 


5.35 ■ 10" 


2 


5.48 ■ 10-' 


0.18 


M2V 


3500 


4.00- 10-' 


3.38 ■ 10- 


2 


5.00- 10-' 


0.20 


M3V 


3500 


3.30- 10-' 


2.74- 10" 


2 


4.50-10-' 


0.21 


M4V 


3500 


2.63- 10-' 


1.64- 10" 


2 


3.49 ■ 10~' 


0.34 



surface brightness {e.g. stellar spots, polar caps for fast rotators) 
or star misalignment with instrumental line of sight. 

Stellar spots, for instance, represent a localized flux de- 
fault (Fig. 1). For a temperature differential of ~ 10 3 K, and 
a spot area of 10% of the star's surface, the flux default is 
50 ph s _1 m -2 fun . However, with an ideal nulling transmission 
on the limb of the star of 10~ 8 , this particular feature will not in- 
terfere with the signal from an exo-Earth. 



3.4. Star Output 

Taking into account the considerations presented in Sects. 3.1 
through 3.3, the option chosen in our model is to generate a 
resolved image of the limb-darkened stellar surface with suffi- 
cient resolution to allow for a correct estimation of the leakage 
(Fig. 2). In order to avoid introducing too many settings in our 



Belu et al.: Exoplanetary scene model 



5 




e 



+ 



Fig. 1. Stellar spot - point source equivalence. When viewed 
through the asymmetric filter that the instrument represents, a 
stellar spot is equivalent to a symmetrically placed virtual planet. 
R is the interferometric coronagraphic transmission, and 9 is an- 
gular sky position parameter. 



□EC (deg) DEC (deg» 




7000 8000 9000 1QO00 62 63 64 6b 



Fig. 2. 209 x 209-pixel Origin output images of a GOV star 
at 10 pc at wavelengths A = 0.5 urn and A = 10 urn. Linear 
grayscales go from to the maximum of brightness in each plot. 
The bandwidth of the channel is 1 [im. 

model (which would prevent inspection of a wide range of pos- 
sible scenarios), the only parameter considered to account for 
non-symmetric stellar contribution is the pointing error, i.e. a 
possible angular offset between the line of sight and the position 
of the star. 

The images in Fig. 2, and a majority of 
the following, are outputs of the widely used fv 
(http://heasarc.nasa.gov/ftools/fv) Fits viewer 
and editor program, which is proposed as an interface for the 
ORIGIN software. This was preferred to a dedicated interface 
development, in order to use tools already familiar to the 
astrophysical community as much as possible. 

4. Exoplanets 

The exoplanets will be unresolved at the instrument's resolution. 
The specific intensity of the j-th planet is therefore: 

Ij(s, A, t) = [F eJ (A, t) + F rJ (A, t)] d(s - Sj(t)) , (19) 

where 6 is Dirac's distribution, Sj{t) is the angular position of 
the planet at date f, F e is the specific flux intrinsically emitted 
by the planet and F r is the stellar flux reflected by the planet. 
Absorption of the planetary flux by exozodiacal dust, even in 
the case of edge-on systems, would assume a global dust density 
that would rule out nulling detection, so this is not implemented. 



The specific fluxes emitted and reflected by a given planet are 
discussed in the following section; we have, however, dropped 
the planet index j for the sake of simplicity. 

4.1. Planet Emission 

The diversity of exoplanets is expected to be considerable, espe- 
cially in the case of terrestrial planets (Gaidos & Selsis 2007). 
The composition of a given terrestrial exoplanet's atmosphere, 
and thus its spectrum, will depend on many parameters. Among 
these parameters are the stellar type and the detailed chemical 
composition of its parent star, its orbital distance, its mass, its 
accretion history, the relative abundance of accreted solids (sil- 
icate, metal, ice) and volatiles (including water), the age of the 
system, and possibly, the existence of an extensive biosphere. 
Simulation of planetary atmospheres is a rather young field, and 
the encompassed models represent a level of complexity beyond 
the scope of this work. The best approach for an exoplanetary 
systems FOV model is to be able to use external spectra, pro- 
vided by teams working on the computation of realistic synthetic 
spectra. A collection of identified object spectra is already avail- 
able: phase-dependent Jupiter-like planets at different orbital dis- 
tances (Barman et al. 2005), Earth, Mars or Venus-like planets 
and derived terrestrial planets (Selsis 2000; Schindler & Kasting 
2000; Selsis et al. 2002; Des Marais et al. 2002; Tinetti et al. 
2006, 2005), Earth-replicas orbiting around G, F, and K stars 
(Selsis 2000; Segura et al. 2003), M stars (Segura et al. 2005) 
the Earth throughout its history (Kaltenegger et al. 2007), Earth- 
like planets across the habitable zone (Paillet 2006), and "ocean- 
planets" (Leger et al. 2004). 

Alternatively, it is possible to implement a simple black body 
emission spectral energy distribution. In our model, its effective 
temperature can be either calculated upon a radiative equilibrium 
with the star, or fixed by the user. The effective temperature of 
the planet T p derived from radiative equilibrium with the star is: 

r p = (l-A p ) 1/4 J^T*, (20) 

where A p is the Bond albedo, and the distance r from the planet 
to the star depends on the observation date t and is computed by 
solving the orbital equations. If a simple black body emission is 
used, the emitted flux is 

F e (A,t) = n{^\ B A (T V ), (21) 
where B A is the Planck function. 

4.2. Reflected Flux 

Depending on their albedo and on their phase with respect to the 
observer, the planets partially reflect the flux of their hosting star. 
The specific flux received by a planet from its star, at a distance 
r from the star is: 

F A {r) = 2n(^ I\]Q A , (22) 
where the factor Q A depends on the star limb-darkening law: 

2,= f DMvdfi. (23) 
Jo 
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Fig. 3. Full night side (solid) and full day side (dashed) emis- 
sions. Night: emission model of a l.lRj up iter planet of albedo 
0.1, orbiting at 0.02 AU from parent star from Barman et al. 
(2005). Day: added full-disk reflection of Kurucz model 120 
(GOV) stellar flux, from Van Hamme (1993). All seen at lOpc. 
Channel width: 1 um. 

The factor Q A can be computed for the considered limb- 
darkening laws: 

' 1/2 (black body) 

(3 — x A )/6 (linear law) 

(9 - 3 X A + 2 y A )l 1 8 (logarithmic law) 

(15 — 5^ — 3 y A )/30 (square-root law) 

The flux reflected by the planet then is: 

F r (A,O = KW,0 f^j] F A (r). 



(24) 



(25) 



F A (r) is the specific flux received by the planet from its hosting 
star as given by Eq. (22) and d is its distance from the observer. 

£) is the phase-dependent reflectivity of the planet, as seen 
at an angle £ from the star. This general expression of the re- 
flectivity enables our model (which is open) to be updated with 
expressions taking into account specular reflection on a planet's 
surface, as well as Mie and Rayleigh diffusion, which is be- 
yond the scope of our present work. For now, we have used a 
Lambertian model Ref(A,^) = A A (p(^). A A is a constant and <p(£) 
is the fraction of the planetary disk surface illuminated by the 
star as seen by the observer (phase): 



m = \ [l + sin(7r-£)] = i 



[1 + sin(v(f) + <o(t)) sin z] (26) 



where i is the inclination of the orbit, v is the true anomaly, and 
a> is the argument of the periastron of the considered planet at 
the time of observation. 

Figure 3 exemplifies the importance of this phenomenon as 
simulated with our model. Actually, the emission spectrum of 
planets so close to their star is not uniform with their apparent 
phase (Harrington et al. 2006). 




Fig. 4. How the optical depth, for a target at a given ecliptical 
latitude b, increases with the target's L2 -relative ecliptical lon- 
gitude a (S-Sun, L-L2 point(Sun - Earth), A-Antisolar direction, 
/-ecliptical longitude). 



5. Other Sources 

5.1. Local Zodiacal Cloud 

In order to compute the local zodiacal dust emission, we inte- 
grate along the line of sight (Fig. 4), from the observer's loca- 
tion towards the target, through a 3-dimensional sampling of the 
Kelsall model (Kelsall et al. 1998). Our approach is thus similar 
to that of Landgraf & Jehn (2001), with some minor differences, 
and from a "background drift", fixed pointing direction perspec- 
tive. Since there is no theoretical outer limit to the exponential 
law of the dust density in this model we integrate to the outer 
limit of the physical cloud; i.e. its collisional origin in the aster- 
oid belt, at SA * 5 AU. 

Given the size of the interferometer's FOV, and with respect 
to the value of this outer limit, the local zodiacal flux contributes 
uniformly to the image. This is equivalent to a global noise level 
in the nulling data, depending only on the target's sky position 
and on the instrument's position at a given observation date. 

Since scattering by dust is negligible compared to the ther- 
mal emission at 10 um, it can easily be seen that it is best, at 
a given mission time, to observe targets in the antisolar ecliptic 
meridian, because that is where the optical depth of the cloud 
is minimal. However, during spectroscopic observations, which 
can be 4 months long, the optical depth along the line of sight 
will vary. 

Let s* be the direction of the target star (the time parameter 
t is omitted for the sake of readability): 



/z(S*,^) 



J-L2T 
' 



J z (r L2 + us i ,,A)du 



(27) 



where the position r = Il2 + u s* is relative to the Sun, Yu_ is the 
observer position, u is the distance from the observer towards 
the source along the line of sight, and where J z is the emitted 
specific intensity per unit of depth, i.e. the power emitted by the 
dust at position r, per unit of volume, per unit of solid angle, per 
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Fig. 5. Local zodiacal cloud emission drift at L2 point over an 
evenly spread interval of 6 months, steadily increasing from the 
anti-solar direction (lowest) to solar (highest). The ecliptical lat- 
itude of the target is 10°, channel width is 1 urn. FOV radius is 
076. 



unit of spectral bandwidth, given by the model of Kelsall et al. 
(1998). 

The direction of observation s* is (RA, DEC) in equatorial 
coordinates and, to perform the integration, it must be converted 
into ecliptic coordinates (/, b). 

Figure 5 shows the drift of the local zodiacal emission back- 
ground, for a given target, from the L2 point, over a evenly 
spread period of 6 months. At around 14 urn the zodiacal drift 
over 3 months 2 is 200%. Thus the global noise level drift is 
as 14 % in that particular channel, over a duration not unlikely 
for a spectroscopic observation. 

5.2. The Exozodiacal Cloud 

In our model, we assume that the exozodiacal cloud dust has 
similar properties to the solar zodiacal cloud (Sect. 5.1), hence 
we simulate it by scaling the model of Kelsall et al. (1998). 
All parameters are free, enabling for instance to produce large 
clumps of dust useful for robustness testing of planet signal ex- 
traction algorithms. 

5.3. Galactic Interstellar Medium 

Emission levels. Galactic interstellar medium (ISM) IR emis- 
sion can reach hundreds of MJy sr~' in the infrared (Schlegel 
et al. 1998). Over the typical 076 field of view of interferometer, 
this equates to 2500 phs~' irT 2 urrT 1 at 10 urn. We carried out 
simulations on ESA code (den Hartog 2006), showing that the 
additional noise level provided by a such a background doubles 
the detection time of an exo-Earth. 



2 Collectors can point no further than ~ 80° from the antisolar direc- 
tion to enable solar shielding 
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Fig. 6. Histogram of galactic interstellar medium emission at 
12.5 (im, for 609 targets form Kaltenegger et al. (2006) 



We have therefore extracted from the Improved Reprocessing 
of the IRAS Survey (IRIS) maps (Miville-Deschenes & Lagache 
2005) 12.5 and 25 urn background emission levels for the target 
list of Kaltenegger et al. (2006). Figure 6 shows the histogram 
for the emission at 12.5 urn. Background emissions (again, over 
the FOV) do not exceed 900 ph s _I irT 2 uirT 1 , with a mean value 
of ~ 41 and a median value of ~ 32. 

It is worth noting that the resolution of the IRIS maps is 3f8. 
How do we know that in such a pixel, we do not include the 
flux of unresolved sources (galactic stars, background galaxies, 
etc.), thus overestimating the flux of the ISM that would actu- 
ally be intercepted by the interferometer's FOV? Typical star 
count values for magnitudes greater that 20 (K band) range from 
3 [log Afdeg~ 2 0.5mag _1 ] for the galactic disk (not in the galac- 
tic plane) according to Girardi et al. (2005) to 4.5 for the bulge 
(Rodgers et al. 1986); this is equivalent to 25-50 stars per 0.5 
mag per IRIS pixel. We can therefore keep in mind that the val- 
ues of the ISM we have extracted may be conservative. However, 
given the above-mentioned statistics on current candidate lists, 
we have not proceeded to implement this correction yet. 

Spectral variability. As can be derived from Verstraete et al. 
(2001, Fig. 1 therein), variations of the continuum emission 
of the ISM, around 18 urn and integrated over the 076 FOV 
of the instrument, can reach a photon noise contribution of 
50 ph s _1 irT 2 urn . ISM background calibration is thus required 
for spectroscopy. 

The reddening absorption due to interstellar molecular 
clouds between ourselves and the targets, given their relative 
proximity, is too faint to be a bias. 

Variability over the FOV. Diffuse background ISM structure 
studies have a resolution limit of 10"(Ingalls et al. 2004). As can 
be seen in Fig. 3 of this reference, extrapolation of the a - -3.5 
power law to the FOV spatial frequency (076 <-> 1.6arcsec _1 ) 
leads to an extrapolated power level of 10~ 7 MJy 2 sr _1 , corre- 
sponding to a statistical flux variation amplitude over a 076 FOV 
of 10~ 9 Jy at 24 urn, or 6 x 10~ 4 ph s -1 trT 2 urrT 1 . This is compa- 
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Fig. 7. Integrated field galaxy counts. Plot has similar profile at 
60 and 170 urn, with maximums steadily decreasing with the 
wavelength. 



rable to the flux of an exo-Earth (typically 0.5 phs -1 irT 2 urrT 1 
at 10 urn). The emission of the cold ISM at 10 (im should be even 
lower, so the variation of the ISM emission in the FOV should 
not be visible in the nulling data processing, at exo-Earth detec- 
tion level. Finally, it can be noted that the parallax of the closest 
targets (072) is only a fraction of the FOV, so there should be no 
galactic ISM background drift with the parallax. 

5.4. Faint Background Objects 

Figure 7 displays an integration of a previously established 
galaxy count histogram (Lagache et al. 2004). As mentioned 
above, an exo-Earth at 10 pc emits 0.5 ph s _1 irT 2 uirT 1 at 10 urn. 
We see that, statistically, there should be 0.01 objects brighter 
than that in a 076 FOV. In order to examine galactic stars in the 
FOV, we can investigate further the reference cited in Sect. 5.3 
(Girardi et al. 2005, Fig. 6 therein). Let us consider the star count 
peak in the K band (2 um, longest wavelength available in this 
reference; in the FOV, that represents 2.2 x 10~ 5 objects per 0.5 
mag bandwidth, so it is not even necessary to consider the in- 
tegrated count. We conclude that it is unlikely that one will en- 
counter background stars emitting at the level of an exo-Earth or 
greater in the FOV. Were faint background objects still to occur 
in the FOV, they could be easily ruled out spectroscopically, and 
by monitoring their orbital motion during revisits. 

In addition to the final remark of Sect. 5.3 concerning the 
ISM emission, the numbers above remain negligible if an ex- 
tended FOV (taking into account the parallax) is considered in- 
stead, so we have not accounted for time variability of 4 g (galac- 
tic background drift). 

Because of the above, our model does not specifically imple- 
ment FOV background stars or galaxies. 



6. Discussion 

It is interesting to consider whether one or several of the source 
features identified above was already an "astrophysical limita- 
tion" in achieving the scientific objective of these exo-Earth find- 
ing missions. We recall that FOV source features where screened 
on the intuitive criterion whereby "anything producing a signal 
at the level of an exo-Earth deserves modeling", in the prospect 
of full end-to-end testing simulations. 

Currently, target catalogs for these survey missions have 
been started; for now, these are based on astrobiological inter- 
est and the few informed astro-engineering requirements avail- 
able today, such as the existence of a secondary star too near and 
bright (Kaltenegger et al. 2006) 3 . Surveys of exozodiacal dust 
levels of nearby stars are planned or in progress 4 ; the statistics 
of this unknown parameter alone may impose dramatic revision 
of the scale and cost of missions, even though elements in favor 
of optimism do exist (Beichman et al. 2006). 

In this context, the features of the sources described here 
are, in most cases, an order of magnitude below the exozodiacal 
cloud issue. Beyond the classical limitation of shot noise from 
unsought-for sources, additional requirements on the integration 
time (hence on the cost) may arise from signal processing speci- 
ficities, such as systematic multiple planet signal disentangling, 
regardless of their relative positions and kinematics. This is the 
purpose of our ongoing work. 

7. Summary 

We have defined a FOV physical model of exoplanetary system 
scenes, and proposed an associated Fits input and output format. 
The input format is a tentative standard for defining exoplane- 
tary systems. The output format is an input format for exoplanet 
seeking instrumental simulators. They are both described in the 
Appendices. 

Submitting a fixed resolution image to an instrumental sim- 
ulator is not practical, given the dynamics in resolution and flux 
between the various typical sources. The output flux is mod- 
eled by a "layering" of various sources (star, planets, dust). Each 
"layer" is of one type among i) a 3-D spectral flux image (for 
resolved sources), ii) a spectral flux and a position (for unre- 
solved sources), iii) a FOV-uniform spectral contribution. For 
each source we have examined the detectable features that need 
to be modeled. 

Depending on the spectral type of the parent star and of the 
wavelength, omitting to model limb-darkening is shown to in- 
duce a bias in the estimation of the leakage noise of up to ~ 35%. 
Stellar spots produce a fainter signal than that of an exo-Earth. 
Local zodiacal drift is found to be smaller than 14 % in any given 
spectral channel. Galactic ISM background for a current list of 
~ 600 targets peaks at 900 phs 1 m" 2 urn -1 , with a mean value 
of ~ 41 and a median value of ~ 32 for a typical FOV of 076. 
Simulations show that a doubling of the detection time of an 
exo-Earth is induced by a 2500 phs 1 irT 2 urrT 1 galactic emis- 
sion background (compared to no background emission, this is 
the order of magnitude of the highest emissions). Finally, back- 
ground galactic stars and distant galaxies as bright as, or brighter 



3 Also, recent unpublished work by Brown et al. 
http : // sco . stsci . edu/tpf_tldb/downloads/ 
TPF_SWG_presentation_f eb24_04 . pdf 

4 Traub and Kuchner, Shared Risk Keck Nuller observations 
http : //planetquest . jpl . nasa . gov/Navigator/ 
keck_sharedrisk. cfm#traub 
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than an exo-Earth are unlikely in the FOV, even extended by par- 
allax drift. 

The model specifications have been embedded by Starlab 
& Thales Alenia Space (formerly Alcatel Alenia Space, see ac- 
knowledgments) into a Java simulator called Origin, soon to be 
open-source, using the input/output definition standards detailed 
in the Appendices. We are in the process of building upon this 
work to obtain an end-to-end simulation approach. 
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Appendix A: Input format 

We have developed a Fits (Hanisch et al. 2001) standard spec- 
ifying the input parameters for modeling an exoplanetary sys- 
tem. This format takes advantage of the building block structure 
of Fits files. This building block approach also enables modu- 
lar storage of stereotypes, as well as the possibility of linking to 
exterior databases in the future. 

A. 1 . Overview and Primary HDU 

All the data necessary to define a FOV (that can be read by 
flux calculators such as the Origin software - Sect. 7) are 
stored in a Fits file. Databases, such as chromatic specific in- 
tensities of bodies, are stored in Fits binary tables (i.e., with 
XTENSION='BINTABLE'), whereas scalar parameters are stored 
in the headers of the extensions. The primary Header Data Unit 
(HDU) of an Origin input data file is informational only and con- 
tains no data 5 . An example of such a primary HDU is provided 
in Table A.l. 

An Origin input data file provides the following HDUs: 

- general scenario parameters (spectral channels, observation 
epoch, duration of the observation and number of snapshots 
to be produced during that time frame, FOV resolution and 
sizes), 

- star parameters, 

- local and exozodiacal cloud parameters, 

- planet(s) parameters. 

The order of HDUs is irrelevant, the file contains at most one 
HDU of each type (except for planets). HDUs are identified by 
their names (value of EXTNAME Fits keyword). Table A. 2 lists 
the different Fits extensions (actually Fits binary tables) used 
to implement the Origin input file format. The revision num- 
ber of the Origin input file format described in this document is 
EXTVER=1 and is indicated by the value of Fits keyword EXTVER 
in each HDU extension. The various extensions used in Origin 
input data files are detailed in the subsequent subsections. 

A.2. Scenario Table 

The parameters defining the scenario of an exoplanetary system 
observation are stored into a Fits extension named ' SCENARIO ' . 
The corresponding binary table contains a first column with 
the central wavelength of the channels, and a second with their 
width. A number of scalar parameters are also provided in the 
header part of this HDU: 

5 According to Fits standard, primary header can only contain image 
data, not binary tables. 
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Table A.l. Example of the header part of a primary HDU in an Origin input file. 



SIMPLE = T / true FITS file (\protect\vrule width8pt\protect\href {http://fits.gsfc.nasa.gov/Hhttp: //fits. gsfc.nasa.gov/}) 

BITPIX = 8 / 8-bit twos complement binary unsigned integer 

NAXIS = 6 / this HDU contains no data 

EXTEND = T / this file may contain FITS extensions 

BLANK 

COMMENT This is an input FITS file for ORIGIN software for exoplanetary 

COMMENT system sky energy distribution modeling. 

BLANK 

HISTORY Created by SOMEBODY on SOMEDATE. 
END 



Table A.2. Description of Fits extensions of Origin input file 
format. All these extensions are saved into a Fits binary ta- 
ble identified by its name, which is the value of the keyword 
EXTNAME. Column Number indicates the number of extensions 
of a given type allowed in the file. 



Table A.3. List of Fits keywords used to define scalar parame- 
ters of the ' SCENARIO ' extension. 



Extension Name 


Number 


Description 


'SCENARIO' 


1 


global parameters 


'STAR' 


1 


star parameters 


'EXO-ZODI' 


1 


exozodiacal cloud parameters 


'LOCL-ZODI' 


1 


local zodiacal cloud parameters 


'UNRESOLVED' 


any 


parameters for planets 




or other point- 






like objects 





- The value of L2-FLAG specifies whether the observer's posi- 
tion is at the L2 point (see Sect. 2.2). Otherwise, Earth posi- 
tion is assumed. 

- The value of DW- EPOCH is the Julian date of the beginning of 
the operational phase of the instrument, as an epoch that later 
mission events will be relative to (DW comes from Darwin, 
the mission for which this standard was initially developed). 

- The value of OBS-DATE is the time (in fractional JD) from 
DW-EPOCH when the current observation starts. Origin uses 
compact computer algorithms by Fliegel & van Flandern 
(1968) for converting between Julian days and Gregorian 
calendar dates. 

- The value of OB S- STEP is the duration between two suc- 
cessive snapshot outputs of the scene, for following orbital 
motion of planets. 

- STR-RES is the resolution at which the flux calculator (as the 
Origin software) generates a chromatic sky energy distribu- 
tion (an image cube) of the star. 

- IMG-RES and IMG-FOV are the resolution and FOV, respec- 
tively, of the image (cube) of the exozodiacal dust. Also, 
IMG-FOV is used to compute the FOVs solid angle for uni- 
form contribution calculation (see Sect. B.4). 

- CLOUD-DZ is the integration step along the line of sight 
through 3-D dust distributions. 

These parameters are listed in Table A. 3. Table A.4 shows 
a typical SCENARIO header of an Origin input file. Note that, 
in this header, the value of NAXIS2 is also the number of effec- 
tive spectral channels for which output will be generated. The 
columns of the binary table in a SCENARIO extension are listed 
in Table A. 5. 



A.3. Star Parameters 

The star model parameters are stored in a 'STAR' extension. 
Table A. 6 lists the keywords of the star model. 



Keyword 


Description 


Units 


L2-FLAG 


observer at L2 point 




DW-EPOCH 


date of beginning of mission 


JD 


OBS-DATE 


observation date (from DW-EPOCH) 


JD 


OBS-STEP 


time between 2 snapshot outputs 


his 


OBS-NB 


number of snapshot outputs 




STR-RES 


star resolution 


// 


IMG-RES 


resolution of the zodiacal image 


ft 


IMG-FOV 


FOV 


ft 


CLOUD-DZ 


z integration step through cloud 


AU 



Table A.5. Description of Fits binary table for the scenario ex- 
tension of the Origin input file format. 



Column 


Description 


Units 


SPCH-CWL 


central wavelength of channels 


m 


SPCH-WDT 


bandwidths of channels 


m 



Origin provides a database of monochromatic limb- 
darkening parameters and specific intensity for stars of various 
spectral type and luminous class (see Sect. 3.1). This database is 
built from Van Hamme (1993) tables and from a model of the HR 
diagram of existing stars, to establish the relation between spec- 
tral type and luminosity class, and physical parameters such as 
the star's effective temperature, surface gravity and luminosity. 
Tables A. 7 displays, for reference, the mean specific intensity in 
photometric bands, for the spectral types of stars envisioned for 
exoplanet search. 

The star model can be built by choosing one of the items 
from the stellar database provided with the Origin software, or 
by specifying the star parameters (either in the same format as 
the database or, more simply, by a simple black body emis- 
sion model characterized by the star effective temperature, lumi- 
nosity or radius, mass or surface gravity and, optionally, limb- 
darkening parameters). 



A.4. (Exo)zodiacal Cloud 

Exo- and local zodiacal cloud parameters are stored into 
EXTNAME=' EXO-ZODI' and EXTNAME= ' LOCAL -ZODI ' exten- 
sions. Table A. 8 lists the keywords of the (exo)zodiacal model. 
Most parameters are those of the implemented Kelsall et al. 
(1998) model. The binary table of these extensions list solar sys- 
tem's Kelsall et al. (1998) parameters of the three dust bands of 
the zodiacal cloud. 
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Table A.4. Typical SCENARIO header in an Origin input file. 



XTENSI0N= 
BITPIX = 
NAXIS 
NAXIS1 = 
NAXIS2 = 
PCOUNT = 
GCOUNT = 
TFIELDS = 
EXTNAME = 
EXTVER = 
TFORM1 = 
TTYPE1 = 
TUNIT1 = 
TF0RM2 = 
TTYPE2 = 
TUNIT2 = 
L2-FLAG = 
DW-EPOCH= 
OBS-DATE= 
OBS-STEP= 
OBS-NB = 
STR-RES = 
IMG-FOV = 
IMG-RES = 
ORB-FOV = 
ORB-RES = 
CLOUD-DZ= 
END 



/ FITS 3D BINARY TABLE 

8 / Binary data 

2 / Table is a matrix 

16 / Width of table in bytes 

13 / Number of entries in table 

8 / Random parameter count 

1 / Group count 

2 / Number of fields in each row 
/ Table name 

1 / Version number of table 

/ Data type for field 

/ Label for field 

/ Physical units for field 

/ Data type for field 

/ Label for field 

/ Physical units for field 

/ Observer at L2 point? 



2456686 / 
166 / 
24 / 
1 

8E-6 
6.715 
6.61 

2 

6.61 
8.1 / added AU 



Beginning of mission, JD 
Beginning of observation, 



from DW-EPOCH 



Table A.6. Star model keywords. 



symbol 


description 


units 


ICRS- ALP, BET 


ICRS star coordinates 


deg 


SHIFT- ALP, BET 


offset of the star image from the pointing direction of the instrument 


deg 


SP-MODE 


whether black body or provided spectrum should be used 


flag 


LIMB-DRK 


which (uniform, linear, quadratic, etc.) limb darkening law should be used 


flag 


RLT-FLAG 


which among luminosity, star radius or effective temperature should be calculated from the others 


flag 


MGR 


which among mass, surface gravity or radius should be calculated from the others 


flag 


ML 


which from mass or luminosity should be calculated from the other 


flag 


DISTANCE 


distance from observer 


pc 


EFF-TEMP 


star effective temperature 


°K 


MASS 


mass of the star 


M B 


LUMINOSI 


star luminosity 


Ad 


RADIUS 


star radius 


l< 


LOG-GRAV 


log of star surface gravity 


go 



A. 5. Unresolved Sources 

Extension(s) EXTNAME= ' UNRESOLVED ' store the parameters for 
unresolved sources such as planets. The format allows the speci- 
fication of all their orbital parameters, but can also be used to 
account for fixed background point-like sources such as field 
stars or distant galaxies. Table A. 9 lists the keywords of the un- 
resolved sources model. In addition to these, Kelsall et al. (1998) 
parameters of the circumsolar dust ring and blob can be speci- 
fied. Origin uses the Kelsall et al. (1998) model equations, ex- 
cept the radius of the ring is the semi-major axis of the planet. 
Table A. 10 lists the columns of the binary table of the unresolved 
source extension: this is the optionally provided planetary spec- 
trum. 

Appendix B: Output Format 

The output format of the Origin software is also proposed as an 
input standard for instrumental simulators. It consists of two Fits 



Table A.10. Description of Fits binary table for unresolved ex- 
tension in the Origin input file format. 



Column 


Description 


Units 


LAMBDA 


channel's central wavelength 


m 


SP_FLUX 


specific intensity 


u -1 -2 -1 -1 

ph s m _ sr urn 



files (or two series of files, if several observations at different 
times are demanded). The ORBIT file simply contains an image 
of the planets' orbits with the current position of each planet, and 
it mainly serves a human visual check purpose. The LAYERED file 
contains the physical information that an instrumental simulator 
needs, and is described in the following. 

We first note that the dynamic in specific intensity and reso- 
lution between the different sources in a planetary system scene 
is such that it is impractical to compute a global image, with a 
given resolution, of the whole FOV, and submit it to an instru- 
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Table A.7. Mean specific intensity in photometric bands (units: ph s 1 m 2 sr 1 |xm '). 



Type 




U 




B 




V 




R 




I 




J 




H 




K 




L 


M 




N 




Q 


FOV 


1.5 


10 26 


2.6 


10 26 


2.4 


10 26 


1.8 


10 26 


1.2 


10 26 


7.1 


10 25 


3.9 


10 25 


1.9 


10 25 


5.5 


10 24 


1.8 


10 24 


2.4 


10 23 


3.9 


10 22 


F2V 


1.3 


10 26 


2.2 


10 26 


2.1 


10 26 


1.6 


10 26 


1.1 


10 26 


6.7 


10 25 


3.7 


10 25 


1.8 


10 25 


5.4 


10 24 


1.8 


10 24 


2.3 


10 23 


3.7 


10 22 


F5V 


8.6 


10 25 


1.5 


10 26 


1.6 


10 26 


1.3 


10 26 


9.6 


10 25 


5.8 


10 25 


3.4 


10 25 


1.7 


10 2S 


5.0 


10 24 


1.6 


10 24 


2.1 


10 23 


3.5 


10 22 


F8V 


6.7 


10 25 


1.2 


10 26 


1.3 


10 26 


1.1 


10 26 


8.7 


10 25 


5.4 


10 25 


3.3 


10 25 


1.6 


10 25 


4.8 


10 24 


1.5 


10 24 


2.0 


10 23 


3.3 


10 22 


GOV 


5.1 


10 25 


9.6 


10 25 


1.1 


10 26 


1.0 


10 26 


7.8 


10 25 


5.0 


10 25 


3.1 


10 25 


1.5 


10 25 


4.6 


10 24 


1.5 


10 24 


2.0 


10 23 


3.2 


10 22 


G2V 


3.7 


10 25 


7.5 


10 25 


9.3 


10 25 


8.6 


10 25 


6.9 


10 25 


4.6 


10 25 


3.0 


10 25 


1.4 


10 25 


4.4 


10 24 


1.4 


10 24 


1.9 


10 23 


3.1 


10 22 


G5V 


3.7 


10 25 


7.5 


10 25 


9.3 


10 25 


8.6 


10 25 


6.9 


10 25 


4.6 


10 25 


3.0 


10 25 


1.4 


10 25 


4.4 


10 24 


1.4 


10 24 


1.9 


10 23 


3.1 


10 22 


G8V 


2.5 


10 25 


5.7 


10 25 


7.6 


10 25 


7.4 


10 25 


6.0 


10 25 


4.2 


10 25 


2.8 


10 25 


1.4 


10 25 


4.2 


10 24 


1.3 


10 24 


1.8 


10 23 


2.9 


10 22 


KOV 


1.7 


10 25 


4.3 


10 25 


6.0 


10 25 


6.1 


10 25 


5.2 


10 25 


3.8 


10 25 


2.6 


10 25 


1.3 


10 25 


4.0 


10 24 


1.2 


10 24 


1.7 


10 23 


2.8 


10 22 


K1V 


1.0 


10 25 


3.1 


10 25 


4.7 


10 25 


5.0 


10 25 


4.4 


10 25 


3.4 


10 25 


2.5 


10 25 


1.2 


10 25 


3.8 


10 24 


1.1 


10 24 


1.6 


10 23 


2.7 


10 22 


K2V 


1.0 


10 25 


3.1 


10 25 


4.7 


10 25 


5.0 


10 25 


4.4 


10 25 


3.4 


10 25 


2.5 


10 25 


1.2 


10 25 


3.8 


10 24 


1.1 


10 24 


1.6 


10 23 


2.7 


10 22 


K3V 


6.0 


10 24 


2.1 


10 25 


3.4 


10 25 


4.0 


10 25 


3.7 


10 25 


3.0 


10 25 


2.3 


10 25 


1.2 


10 25 


3.6 


10 24 


1.1 


10 24 


1.5 


10 23 


2.5 


10 22 


K4V 


3.2 


10 24 


1.4 


10 25 


2.4 


10 25 


3.1 


10 25 


3.0 


10 25 


2.6 


10 25 


2.2 


10 25 


1.1 


10 25 


3.4 


10 24 


9.9 


10 23 


1.5 


10 23 


2.4 


10 22 


K5V 


1.6 


10 24 


8.1 


10 24 


1.6 


10 25 


2.2 


10 25 


2.4 


10 25 


2.2 


10 25 


2.0 


10 25 


1.0 


10 2S 


3.2 


10 24 


9.2 


10 23 


1.4 


10 23 


2.3 


10 22 


K7V 


7.6 


10 23 


4.4 


10 24 


9.9 


10 24 


1.5 


10 25 


1.8 


10 25 


1.8 


10 25 


1.7 


10 25 


8.8 


10 24 


2.9 


10 24 


8.6 


10 23 


1.3 


10 23 


2.2 


10 22 


MOV 


3.9 


10 23 


2.3 


10 24 


5.8 


10 24 


9.6 


10 24 


1.3 


10 25 


1.4 


10 25 


1.3 


10 25 


7.4 


10 24 


2.6 


10 24 


8.0 


10 23 


1.3 


10 23 


2.1 


10 22 


M1V 


3.9 


10 23 


2.3 


10 24 


5.8 


10 24 


9.6 


10 24 


1.3 


10 25 


1.4 


10 25 


1.3 


10 25 


7.4 


10 24 


2.6 


10 24 


8.0 


10 23 


1.3 


10 23 


2.1 


10 22 


M2V 


1.8 


10 23 


1.2 


10 24 


3.2 


10 24 


5.5 


10 24 


9.4 


10 24 


1.1 


10 25 


1.1 


10 25 


6.2 


10 24 


2.3 


10 24 


7.3 


10 23 


1.2 


10 23 


2.0 


10 22 


M3V 


1.8 


10 23 


1.2 


10 24 


3.2 


10 24 


5.5 


10 24 


9.4 


10 24 


1.1 


10 25 


1.1 


10 25 


6.2 


10 24 


2.3 


10 24 


7.3 


10 23 


1.2 


10 23 


2.0 


10 22 


M4V 


1.7 


10 23 


1.1 


10 24 


3.2 


10 24 


5.7 


10 24 


9.5 


10 24 


1.1 


10 25 


1.0 


10 25 


6.1 


10 24 


2.3 


10 24 


7.4 


10 23 


1.2 


10 23 


2.0 


10 22 



Table A.8. (Exo)zodiacal model keywords. 



symbol description units 



TYPE 


whether this extension is a local or exo- zodiacal cloud 


flag 


ZODI 


mean density of dust 


zodi 


DUST-ST 


dust sublimation temperature 


°K 


DUST-RT 


dust reference temperature 


°K 


ETL 


exponent temperature law coefficient 


unitless 


CLOUD-0 


cloud outer radius 


AU 


SC-* 


smooth cloud parameters of the Kelsall model, not all listed here 





Table A.9. Unresolved sources (planets) model keywords. 



symbol 


description 


units 


ORB-APER 


which, from semi-major axis or period is computed from the other 


flag 


P-RADIUS 


planet radius 


R e 


P-TEMP 


planet temperature 


°K 


P-ALBEDO 


planet albedo 


unitless 


ORB-TO 


epoch of periastron passage 


Julian date 


ORB-OMA 


position angle of the ascending node 


deg 


ORB-OMP 


argument of periastron 


deg 


ORB -I 


inclination of orbit 


deg 


ORB-E 


eccentricity 


unitless 


ORB-A 


semi-major axis 


AU 


ORB-PER 


orbital period 


days 


SP-FLAG 


whether to use specific intensities provided in the extensions binary table, or black body emission 


flag 



mental simulator. Following the input standard, the output is also 
layered, as a Fits file containing three types of descriptions of in- 
coming fluxes: resolved sources, unresolved sources, and FOV- 
uniform contributions, that can all be used by an instrumental 
simulator. 



B.1. Primary HDU 

Unlike the input format, the primary HDU restates some of the 
data from the input format. Table B.l gives an example of such 
a header. 



B.2. Output Spectral Channels 



These files, being generated by the Origin software, are quite 
explicit, so this section is considerably more straightforward 
than the previous. 



In the output file, the EXTNAME= ' SPCHANNELS ' extension is a 
binary table which gives the effective central wavelengths and 
spectral bandwidths of the simulated model. This information is 
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SIMPLE = T / Java FITS : Mon May 22 87:38:11 CEST 2866 

BITPIX = 8 

NAXIS = 8 / Dimensionality 

EXTEND = T / file contains FITS extensions 

COMMENT = 'ORIGIN LAYERED OUTPUT FILE' / 

ORIGIN = 'ORIGIN vl.8' / Name and version of software 

AUTHOR = 'Author ' / 

DATE = 'Mon May 22 87:38:18 CEST 2886' / File creation date 

CTRLFILE= 'C:\ORIGIN\Livraison\Full Scenarii\Gl_876. fits' / Name and path of ORI 

DW-EPOCH= 2456888.8 / Darwin epoch [JD] 

OBS-DATE= 188.8 / Observation date since epoch [days] 

IMG-NBR = 6 / Number of images generated during simulation 

IMG-RANK= 1 / Rank of current image simulation 

END 



Table B.3. Description of Fits binary table columns for constant that the columns of this table are exactly the same as those of 
layer extension in the Origin output file format. the CONSTANT LAYER extension (Sect. B.4 and Table B.3). 



Column 


Description 


Units 


SPCHCWL 


channel's central wavelength 


m 


SPCHWDT 


channel's width 


m 


SPCHFLUX 


flux 


ph s -1 nT 2 



similar to that specified in the input format (see Sect. A. 2), and 
is not further described here. 



B.3. Star or Exo-zodiacal Output Maps 

In the output file, the EXTNAME= ' RESOLVED OBJECT LAYER' 

extensions are Fits images containing 3-D maps (right ascen- 
sion, declination and wavelength) of the star's limb darkened 
photosphere, or of the exozodiacal dust cloud. Table B.2 shows 
an example of the header of this extension. The number of pixels 
in the image is always odd, with the reference pixel marking the 
center of the image. For an image example, refer to Fig. 2. 



B.4. Local Zodiacal and Galactic Background 

The background emission of the local zodiacal light, which is 
uniform over the field of view, is provided by the output for- 
mat in phs 1 irT 2 per spectral channel in the binary table ex- 
tension EXTNAME= ' CONSTANT LAYER' (Fig. 5). For that, the 
Origin software considers the FOV solid angle defined in the 
EXOZODI extension in the input format (see Sect. A. 2). The 
header of this extension contains only the binary tables column 
definitions (Table B.3). 

The galactic background can be manually inserted here using 
this same extension format, since it requires no calculation by 
the Origin software. 



B.5. Planets 

The output file format contains as many 
UNRESOLVED OBJECT LAYER binary table extensions as 
there were UNRESOLVED planet definition extensions in the 
corresponding input file. Table B.4 shows the header of such 
a planet binary table output extension. It provides the precise 
coordinates in the field of view of the unresolved source 
(fields RA-STR and DEC-STR). The binary table itself contains 
the calculated flux for the planet, degraded to the resolution 
indicated by the primary extension (Sect. B.l). It can be seen 
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Table B.2. Typical 'RESOLVED OBJECT LAYER' header in an Origin output file. 



XTENSI0N= 

BITPIX = 
NAXIS 

NAXIS1 = 

NAXIS2 = 

NAXIS 3 = 

PCOUNT = 

GCOUNT = 

EXTNAME = 

EXTVER = 

CTYPE1 = 

CRPIX1 = 

CRVAL1 = 

CDELT1 = 

CTYPE2 = 

CRPIX2 = 

CRVAL2 = 

CDELT2 = 
END 



' IMAGE 



/ Java FITS: Wed Mar 15 16:35:36 CET 2666 



/ Dimensionality 



-64 
3 
121 
121 
13 

/ No extra parameters 

1 / One group 

'RESOLVED OBJECT LAYER'/This Layer contains the star image 
1 / 

'RA ' / Right ascension axis 

61 / Reference pixel along axis 1 (starting from 1) 
75.0 / Coordinate of ref. pixel along axis 1 [deg] 
-2.222222222222222E-9/ Pixel step along axis 1 [deg] 
'DEC ' / Declination axis 

61 / Reference pixel along axis 2 (starting from 1) 
155.0 / Coordinate of ref. pixel along axis 2 [deg] 
2.222222222222222E-9 / Pixel step along axis 2 [deg] 



Table B.4. Typical ' UNRESOLVED ' header in an Origin output file. 



XTENSION: 
BITPIX : 
NAXIS 
NAXIS1 : 
NAXIS 2 : 
PCOUNT : 
GCOUNT : 
TFIELDS : 
TF0RM1 : 
TDIM1 
TF0RM2 : 
TDIM2 
TF0RM3 : 
TDIM3 
EXTNAME : 
EXTVER : 
RA 
DEC 

RA-STR : 
DEC-STR : 
TTYPE1 
TUNIT1 
TTYPE2 
TUNIT2 
TTYPE3 
TUNIT3 
END 



/ Java FITS: Wed Mar 15 16:35:43 CET 2666 

8 

2 / Dimensionality 
24 
13 



1 

3 

' ID 

'(1) 
' ID 
'(1) 
' ID 
'(1) 

'UNRESOLVED OBJECT LAYER'/This Layer contains the spectrum of a planet 

1 / 

75.66665875513867 / Object Absolute Right Ascension [deg] 
155.6661242618757 / Object Absolute Declination [deg] 
6.21151849913666165 / Object Right Ascension from Star [Arcsec] 
6.4473427524966594 / Object Declination from Star [Arcsec] 



= 'SPCHCWL 



' SPFLUX 
' count/m*2/s ' 



/ Spectral channel central wavelength 

/ 

/ Spectral channel width 

/ 

/ Object Flux 

/ 



